o 

OS 
Os 

PL, 



Tunneling between two Luttinger liquids with long range interaction 
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The non linear charge transfer through a tunnel junction between two Luttinger systems is 
studied for repulsive, finite range interaction between electrons on the same, Vii, and on different, 
Via, sides of the junction. Features of the Coulomb blockade effect are observed if V12 = 0. We 
predict a novel interaction induced enhancement of the current if V12 > 0. When V12 = V11, 
the current is suppressed at small bias, but the "charging energy", obtained from the asymptotic 
behavior at high bias voltage, vanishes. 

PACS numbers: 73.40.Gk, 72.10.Bg, 72.10.-d 
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Due to the repulsive interaction between the electrons, 
tunneling through mesoscopic tunnel junctions is sup- 
pressed |I],|| for voltages U < Uc = e/2C, and tem- 
peratures T < Tc = Ec/ks p (fcs Boltzmann con- 
stant, e elementary charge, C capacitance). The quan- 
tity Ec = eUc is called charging energy, the suppression 
of the transport Coulomb blockade [l]] . 

In the semi-phenomenological theory |4|] the tunnel 
junction is replaced by a capacitance and a tunnel re- 
sistance. An "external" impedance Z{uf) represents the 
coupling of the tunneling particles to a reservoir of dis- 
sipative degrees of freedom. When Z(0) = R = 0, 
the current-voltage characteristic I(U) is linear. For 
R =/= 0, the current is depleted for U <C Uc- When 
Uc -^iU^ooIrxU — Uc- The shift Uc is an impor- 
tant feature of the Coulomb blockade phenomenon for 
R — > 00 which is also found in experiments H. 

In the quantum mechanical theory, the Hamiltonian 
of the system, even using the Luttinger approximation 
M, becomes non linear. For zero range repulsive interac- 
tion between the electrons, it was shown by using renor- 
malization group arguments raj?! and by conformal field 
theory [R[ that an infinitesimally small scattering bar- 
rier in a one dimensional (ID) Luttinger liquid becomes 
completely insulating at zero temperature. The current 
is suppressed at small voltage, I(U) ex [Z 2 /^ 1 (interac- 
tion parameter g < 1 for repulsive interaction). This 
result was also found to the order A 2 for the model of 
two Luttinger systems with N (3> 1) transport channels 
connected by a tunnel junction described by a transmis- 
sion probability A [|9| . An effective interaction parameter 
g(N) was identified which tends to unity for N — > 00. 

In the present paper, we study a tunnel junction con- 
necting two ID Luttinger systems with finite range in- 
teraction. The correlations between the particles located 
at different sides of the junction are taken into account. 
The current is calculated in terms of the density-density 
correlation function. We recover the above mentioned 
depletion of the current at small voltages, independent 
of the relative strength of the interaction between elec- 
trons on the same, V11 and on different sides, V12, of 



the junction. However, the latter turns out to determine 
the behavior at higher bias voltages: for V 12 = the 
"Coulomb blockade" is found. For V12 > 0, we find that 
the depletion of the current at high voltage is reduced. 
When V12 — V11, the current approches asymptotically 
the non interacting, linear limit. As a consequence, a 
charging energy can only be obtained for an interaction 
potential with a finite, non zero range, when V12 < V\\. 
For V12 > Vii, though the interaction is assumed to be 
repulsive, the current is larger than without interaction. 
We consider the Hamiltonian H = Hq + H t + Hjj. The 
unperturbed part, Hq 7 describes two separate spinless 
Luttinger systems, 1 and 2. They are assumed to extend 
from — L to and from to L (L — > 00), respectively. 
The tunnel junction (at x = 0) is represented by H tl and 
Hu is the energy contributed by the external voltage. 
The interaction energy is 



■fftet = 2 / dxdy p( x )p(y) 

xV(\x - y\) [\xi&{xy) + Xi 2 @(-xy)] 



(1) 



The interaction potential V(\x — y\), is assumed to in- 
troduce a length scale, say a" 1 . For a screened poten- 
tial, V(x) oc exp (— |x|/£), oT x — £. For a Coulomb 
potential in ID, V(x) oc (x 2 + d 2 )" 1 / 2 , a" 1 = d @. 
The (real and positive) parameters An and A12 are in- 
troduced in order to vary the strengths of Vn = XnV 
and V12 = Ai2^ separately. The density operator is 
p(x) = p^(x)9(-x) + p^(x)e(x) (@(x) Heavyside 
function) . 

A crucial point is that the boundary conditions [[flj are 
such that the original Fermion fields vanish at x = 0, ±L. 
This implies that the corresponding left and right mov- 
ing parts are not independent but ^^(x) = —^ I 



W'i-x), 



rU), 



- iT/0') 



^r( x + 2£) = tyjl(x), and either one of the two alone 



suffices to describe the system 



Then, neglecting 2/cf- 
scattering, p u >(x) = p^ (x) + p^ '(— x), j = 1,2, such 
that the Fourier transform of the interaction Hamilto- 
nian contains terms that are non diagonal in the wave 



numbers. We will see below that these lead to consid- 
erable complications in the calculation of the non linear 
current-voltage relation and the charging energy. 

The tunneling Hamiltonian is given by H t = H^ + 



,( 2 )t, 



r(l)/ 



Hf = £A(*£ ;, (0)*}j ; (0) + h.c.). The energy associated 
with the external voltage is Hu = —ef_ L dxU(x)p(x). 
We assume that the bias voltage U(x) drops accross the 
tunnel contact only ]L2[ . 

The average of the current operator / = ie[Hj~ — H t f] 
is calculated from the backward and forward scattering 
rates 7* = f™ x dt(Hf(t)H? (0)). The average can be 
performed by using the phase fields ^^(x,t) that are 
related to the Fermion fields as usual ]ll| , 

- A 2 POO 

I(U) = e [7+ - 7-] = ¥=- / dt sin (eUt) e~ w ^ (2) 
^ J -00 



where 



W(t) = — <tf*(t)«5*(0) + (<5$(0)) 2 ) = / dw^^ x 



J(w) 



(1 — cos (tot)) coth ( — ) + i sin (cut) 



, (3) 



with 6&(t) = $ (2) (0, i) - $ (1) (0,£) and the spectral cor- 
relation function J(lo). This shows that the non linear 
transport characteristic of the tunnel barrier is related to 
local current fluctuations represented by the phase fields 
$W, i.e. the dynamic and not the static properties of the 
system. This will be seen below in more detail, when we 
calculate the charging energy in terms of the interaction. 
Generally, the spectral function can be decomposed 



J(oS) = Jn(w) + J22M + Ji2(w) + J2i(w). 



(4) 



It is determined by the excitation spectrum of the 
Bosonic bulk modes, w{q) — vp\q\/g(q) with g~ 2 (q) = 1 + 
\nV(q)/irvF. For small wave numbers q, u>(q) = vp\q\/g 
such that J(lu — > 0) is given by the enhancement of the 
Fermi velocity vf with the interaction parameter, g = 
5(0), that contains only An, J(u) — 2uj/g + 0(uj 3 /a 2 ). 
The first two terms in eq. (|J) represent the contribu- 
tions of the separate left and right Luttinger systems, and 
are given by the expectation values ($(•?) (0, t)^) (0, 0)). 
That they dominate the current at small voltages inde- 
pendent of a and A12, is reflected in the behavior of J for 
small lo. 

The off-diagonal terms in eq. (|j) are due to the average 
(^^(Oji^^OjO)). They are typical vertex contribu- 
tions to the average {H t (t)H t (0)}, and describe correla- 
tions between the charge fluctuations on both sides of 
the junction. They influence in a characteristic way the 
spectral function at high frequencies, 



J{lo — > 00) = 2uj 



i- — (\ii(u>v u y + 2\ l2 v u ) 

Ztt 



where V u = V(q — uj/vp). 

Since the interaction range is assumed to be finite in 
position space, the second term on the right hand side 
of eq. (pf) vanishes for to — > 00 such that only the con- 
tribution of the free, non interacting electrons survives, 
namely J(u>) = 2lo. The presence of the term oc A12 on 
the right hand side reflects the above mentioned fluctu- 
ation correlations. We will see below in detail that the 
latter counteract the suppression of the current induced 
by the coupling to the bulk modes and eventually lead 
to a vanishing of the charging energy if An = A12. If 
A12 > An the current at high bias voltage becomes even 
enhanced above the value without interaction. The be- 
havior of the current-voltage characteristic at zero tem- 
perature is shown schematically in Fig. |l|. 
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FIG. 1. The zero temperature current- voltage characteris- 
tic, schematically, of a tunnel junction (tunneling resistance 
Rt) connecting two semi infinite Luttinger liquids for various 
ratios of the interaction strengths between electrons on differ- 
ent, A12, and the same, An, sides of the junction. Full line: 
A12/A11 = 0; short dashes: A12/A11 < 1; dots: A12/A11 > 1; 
long dashes: no interaction. 

The connection with the earlier, semi classical results 
8 is established by rewriting I(U) as 



1 _ p-0eU roo 

I ( U) = —^-l ,,E 



dE'f(E)[l-f(E')} 



xP(E + eU~E'), (6) 



(5) 



with the Fermi function f(E) and 

P(E) = — I" dte iEt e-l w W- w ^W. (7) 

2tt J_ 00 

Only in the semi classical approach, the latter function 
plays the role of a probability density for a bulk excita- 
tion of energy E. In the present microscopic model, it 
can become negative. This is a result of the extraction of 
the Fermi factor /(l — /) in the integrand in eq. (0) which 
is artificial, though formally correct, for the present inter- 
acting system. The tunnel resistance R t = 2oj 2 /e 2 A 2 ir 
contains the cutoff frequency lo c which serves to regular- 
ize the Luttinger model. 

The result for zero temperature, 



I{U) 



1 

eRt 



eU 



dE(eU - E)P(E), 



(8) 



can be used to identify the parameters in the semi clas- 
sical approach in terms of the microscopic model. The 
function that corresponds to the impedance in the semi 
classical theory is Z(uj) = J{u)/ui — 2. For small bias, 
eU <C avp, only the behavior of P{E) for small en- 
ergies is important. It is determined by the diagonal 
contributions in eq. (|J) and always positive. The elec- 
trons loose energy via dissipation into the bulk modes. 
We find I(U) » (U/R t )(eU/av F ) 2/9 ~ 2 , which yields the 
dissipative resistance || R/Rr = Z(0) = 2(g _1 — 1) 
(Rk = h/e 2 von Klitzing constant). This shows that the 
elementary excitations of the Luttinger system play here 
the role of the dissipative Bosonic degrees of freedom of 
the " electromagnetic environment" . 

At high bias voltage, eU > avp, 3\i and J21 become 
dominant. They counteract the two dissipative diagonal 
terms Jn and J22. While the latter always supress, the 
former tend to enhance transport. Formally, this can be 
seen from P(E — > 00). We start by calculating the second 
derivative of the current- voltage relation, cq. (j|) I"(U) = 
eP{eU)/R t . For very large U, P(eU) = Z(eU)/eU. If 
we make the (not very restrictive) assumption that for 
large wave numbers V(q) oc q~ 2 , we have itZ(uj — » 00) = 
2K>(Aii/2 — A12) which changes sign when A i2 > An/2. 
When the coupling between electrons left and right of the 
junction exceeds a certain strength, the curvature of the 
current-voltage curve changes from positive to negative. 
This occurs for voltages that correspond to the intrinsic 
length scale a -1 . 

For U 3> avF, only short times contribute to the in- 
tegrand of eq. (0), the current, eq. (g), becomes I(U) — 
R^iU — Ec/e) + 0(C/ _1 ) with the charging energy 

/■OC 

E C = dujZ(u J ) = V{x = 0)(\ 11 -X 12 ). (9) 

Jo 

The derivation of this result is far from trivial. The start- 
ing point is the relation between W(t), eq. (pf), and the 
fluctuation correlation function G kk > (t) 



4tt 



zZ TU G kk'{t)- 



k,k'>0 



(10) 



In W(t), the average (J$(0, t)5$(0, 0)} enters. It can 
be expressed via the dissipation-fluctuation theorem in 
terms of the imaginary part of the Fourier transform of 
F(t). This leads eventually to the identity 



J(u) = uj 2 lmF{uj) 



(11) 



The set of equations of motion for the correlation func- 
tion can be closed for our model. The resulting lin- 



ear integral equation is used to expand G kk i (w) in pow- 
ers of Gf\u)B kk ,{u), where G k 0) = oj 2 /{u 2 - oj 2 ) is 
the unperturbed correlation function. It contains the 
excitation spectrum of the bulk modes. The function 
Bkh' contains the off-diagonal part of the interaction, 
B kk ,{uj) = -(An + \i2)V{k,k')kk' /lu 2 t:L with 



V(k,k') 



dq 



V(q)q 2 



{q 2 - k 2 ){q 2 - k' 2 )' 



(12) 



By performing first the w-integration (cf. eq. (|9j)) and 
then the summations over the wave numbers, one can 
show that in the integral over the imaginary part of F 
only the term linear in B kk i contributes. All higher or- 
der contributions vanish, due to exact sum rules. The 
calculation leads eventually to 



Er 



du> 



k,k' 



Eq 



2(An + A12) 



kk' 



dkdk'V(k,k'). (13) 



The first term in this expression is the result obtained 
without taking into account the off-diagonal interaction 
terms, eq. (|) ftf], E% = 2\ 11 V{x = 0). The second 
can be further evaluated and gives —(An +Ai2)V(ai = 0) 
such that finally the above result eq. (0) is obtained. 

There are several comments which have to be made 
at this point. When the interaction between electrons 
left and right of the tunnel junction is omitted, An = 1, 
A12 = the preliminary result obtained previously |l3] ] 
is recovered apart from a factor 1/2 which is due to dis- 
carding the influence of the boundary conditions on the 
fields. The latter influence is negligible for small frequen- 
cies |yj. For high frequencies, the boundary conditions 
influence quantitatively the result for the charging en- 
ergy, though the qualitative behavior with the strength 
and the range of the interaction is the same as before. 
When the interaction of electrons in the left and the right 
part of the system is switched on, A12 7^ 0, the charg- 
ing energy is reduced as compared to the previous case. 
For a given voltage, the current is increased, due to the 
presence of the additional, (formally) repulsive electron- 
electron term. In principle, when the interaction between 
electrons in the left and the right part of the system dom- 
inates, An < A12, the current increases even to values 
above that of the non-interacting system, I(U) > U/Rf. 

At the first glance, this seems counter-intuitive. How 
can a repulsive interaction term increase the current? 
The puzzle can be solved by noting that this happens 
at high voltages. For sufficiently small bias, in the regime 
dominated by the coupling to the bulk modes, the cur- 
rent remains depleted well below U/Rf. When the bias 
voltage is sufficiently high, net charges with opposit signs 
are introduced left and right of the tunnel junction. This 



is also suggested by the model of the " Landauer dipole" 
at an impurity in the presence of a stationary current 
Q. As a consequence, the net interaction between the 
left and right Luttinger liquids becomes attractive and 
tends to decrease the charge difference - the bias voltage 
- via additional (tunnel) current. This cannot happen 
in the semi classical model since there the quantum pro- 
cesses are introduced into the model "by hand" , and not 
inherently incorporated into the theory. 

Insisting nevertheless on assigning a capacitance to the 
tunnel junction leads to a dependence of the latter on the 
bias voltage pl|: for small U, the capacitance would be 
finite, due to the suppression of the current well below the 
value of the system without interaction. When increasing 
the bias voltage the correlations between the electrons on 
the left and the right of the junction increase the current, 
i.e. decrease the effective charging energy such that the 
"capacitance" would increase to infinity when A i2 — > An. 
This shows its dynamical origin |10[ . 

A similar result is obtained when considering a poten- 
tial barrier in a Luttinger liquid, where the interaction 
between electrons left and right of the barrier is taken 
into account automatically, An = A12. 

Experimentally, it is very difficult to observe Coulomb 
blockade for a single tunnel junction between two metallic 
due to the presence of a large shunt capac- 
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itance 
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By adding more junctions in series, this shunt 
capacitance is reduced due to the small capacitances of 
the other junctions. Our above result offers a microscopic 
interpretation: adding junctions in series reduces the in- 
teraction between electrons in the left and right leads. As 
a consequence, the observed charging energy is increased. 

A direct experimental test could be performed by using 
quantum wires based on semiconductor hetero structures 
and adjusting the electron density such that only one 
subband is occupied. In the presence of a gate across the 
region of the tunnel barrier, the interaction between the 
leads on the left and the right hand side of the barrier 
could be locally screened such that the charging energy 
should appear. 

In summary, we obtained the non linear current- 
voltage characteristic for a model of two ID quantum 
wires of interacting electrons connected by a tunnel junc- 
tion. The features of the Coulomb blockade phenomenon 
were identified. The charging energy was determined as 
a function of the parameters of the interaction potential. 
We showed that it is crucial in this microscopic model 
that the interaction, besides being of finite range, must 
not be too strong between electrons on different sides of 
the junction in order to produce a charging energy. 
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